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Abstract. Observations with the SZ effect constitute a powerful new tool for inves- 
tigating clusters and constraining cosmological parameters. Of particular interest is 
to investigate how the SZ signal correlates with other cluster properties, such as the 
mass, temperature and X-ray luminosities. In this presentation we quantify these 
relations for clusters found in hydrodynamical simulations of large scale structure 
and investigate their dependence on the effects of radiative cooling and pre-heating. 
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1. Introduction 

Galaxy clusters are the largest gravitationally bound objects in the 
Universe. They have typical masses of 10 14 — 10 15 /i _1 M and contain 
hundreds of galaxies within radius of a few Mpc. The intra-cluster 
medium (ICM) is filled with hot ionized gas, typically at temperatures 
1-15 keV, which produces strong X-ray emission and causes spectral 
distortions in the CMB spectrum via the Sunyaev-Zel'dovich (SZ) 
effect (see Ref. [11]). 

Numerical simulations indicate that the non-baryonic dark matter 
component in clusters, which is the dominant fraction of their mass, is 
remarkably self-similar for systems in approximate state of equilibrium, 
see e.g. Ref. [8]. However, the baryonic gas component does not share 
the same degree of self-similarity. This is more evident from obser- 
vations of the Lx~T relation in clusters, which is much steeper than 
is predicted by simple self-similar scaling laws, specially for low-mass 
systems (see e.g. Refs. [4, 14]). This deviation from self-similarity has 
been interpreted as due to non- gravitational processes, such as radiative 
cooling and heating, that raise the entropy of the gas, see e.g. Ref. [10]. 

The purpose of this study is to investigate the correlation between 
intrinsic properties of clusters found in simulations, and to compare 
them with analytical scaling laws. We use high-resolution hydrody- 
namical simulated clusters to assess the impact of radiative cooling and 
pre-heating on these scaling relations. Here we focus only on scalings 
involving the SZ integrated signal, Y, at redshift zero. The effects on 
the mass-temperature and X-ray scaling relations have already been 
analyzed in detail in Refs. [6, 7] for this same set of simulations. A more 
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detailed analysis of the SZ scaling relations can be found in Ref. [2], 
whose results supersede those presented at this conference. 



2. Scaling relations in Clusters 

The quantity we want to study is the total thermal SZ flux density 
received from a cluster. This is defined as the integral of the thermal 
SZ specific intensity (described by the y Comptonization parameter) 
over the solid angle occupied by the source in the sky. More exactly, 
we want to investigate how the frequency-independent SZ flux, 

Y = jydQ = d A 2 JydA= 1 ^-d A 2 j^T c n c dV, (1) 

correlates with other cluster properties. This is known as the SZ inte- 
grated Y-flux or Y-luminosity. Note that dA is the angular diameter 
distance, dA = d£ld A 2 is the sky projected area of the source and 
the last integral is performed over the volume of the cluster. Since y is 
dimensionless, Y has dimensions of a solid angle. Equation (1) indicates 
that Y is proportional to the mean electron temperature, (T e ), and the 
total cluster mass, M, 

Y oc / gas (T c ) M d A 2 , (2) 

where / gas is the gas mass fraction of the cluster. Note that (T e ) is the 
mass (or particle) weighted temperature, which may differ from the 
emission-weighted temperature measured in X-ray observations. 

The dependence of the X-ray bolometric luminosity on mass and 
temperature follows from its definition, 

L ^ij^ A{T)dV - (3> 

The integral is over the volume of the cluster, fim p ~ 1 x 10~ 24 is the 
mean mass per particle, A(T) is the gas cooling function and p gas is 
the gas density. Assuming the gas is well described by an isothermal 
temperature profile one obtains, 

L x oc / 2 as M A c Pcrit A(T) oc / 2 as M A c p crit T 1 / 2 , (4) 

where the last step results from assuming that the bolometric luminos- 
ity in clusters is dominated by bremsstrahlung emission, A oc T 1 / 2 . 

As suggested by the virial theorem the mass and the temperature 
of a cluster are tightly correlated quantities. At the virial radius these 
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are related by the virial relation, /cbT cx GM/r. The mass enclosed 
within this radius can be written as M = Air r 3 A c p C rit/3, where A c 
is the overdensity contrast between the mean cluster density within r 
and the critical density, /? C rit- Using this expression one obtains 

T oc M 2 / 3 (A c Pcrit ) 1/3 oc ( A-) M 2 / 3 (1 + z) . (5) 

Using Eq. (5) into Eqs. (2) and (4), with (T e ) ~ T, one can de- 
rive general scaling relations for the SZ and X-ray luminosities as a 
function of cosmology, redshift, and the cluster mass (or temperature). 
According to these, clusters are expected to present the following inter- 
dependencies of mass, temperature, X-ray luminosities and integrated 
SZ fluxes at redshift zero: 

TocM 2 / 3 , L x ocM 4 / 3 ocT 2 , Y oc M 5//3 cx T 5 / 2 oc L 5 J- 4 . (6) 

These expressions coincide with the scaling relations predicted by the 
self-similar model [5], which is expected to describe well the correla- 
tion between the above properties if shock heating from gravitational 
collapse is the dominant mechanism driving the gas evolution. 

3. Simulations, Cluster identification and catalogues 

We present results from three simulations of a single ACDM cosmology 
with density parameter, f2 m = 0.35, cosmological constant, Q\ = 0.65, 
Hubble parameter h = 0.71, baryon density, Ob^ 2 = 0.019, shape 
parameter of the CDM power spectrum, T = 0.21, and normalization 
as = 0.9. We use 160 3 particles of each gas and dark matter in a box of 
side 100 /i _1 Mpc, which gives particle masses of m gas = 2.6 x 10 9 /i -1 M Q 
and m-dark = 2.1 x 10 10 /i _1 M respectively. The gravitational softening 
at redshift zero was 25 fo _1 kpc. The simulations were carried out as part 
of the Virgo Consortium program, using a parallel version of the hydra 
N-body/hydrodynamics code (see Refs. [1, 9]). 

The first of the three simulations was performed without any ad- 
ditional heating or cooling mechanisms; we refer to this as the 'non- 
radiative' simulation. The second simulation included a model for ra- 
diative cooling using the method described in Ref. [13], except that 
we adopted the cooling tables of Ref. [12] and a global gas metallicity 
evolving as Z = 0.3(i/to)^o, where Z & is the solar metallicity and t/to 
is the cosmic time in units of its present value. Our third simulation, 
which also includes cooling, was performed in order to study the ef- 
fects of pre- heating the gas at high redshift. We considered an energy 
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injection of 0.1 keV per gas particle at redshift four for this run. Note 
that this is less energy injected per particle than that used in the pre- 
heating run of Ref. [2] (1.5 keV). We used identical initial conditions to 
enable direct comparisons between clusters forming in the three runs. 
For more details see Ref. [3]. 

In simulations clusters are identified with a spherical overdensity 
group finder algorithm, as described in Ref. [6]. The identification 
process consists of several steps. The first is to create a minimal- 
spanning tree of all dark-matter particles whose density exceeds 5 = 
178 x Q~°' 55 (z) times the mean dark-matter density in the box. The 
next step is to split the minimal-spanning tree into clumps of particles 
using a maximum linking length equal to 0.5(5 _1//3 times the mean 
inter-particle separation. Then we grow a sphere around the centre 
of the clump until the enclosed mean density is A times larger then the 
comoving critical density. Because the density parameter of the real 
universe is not known, we have chosen to average the properties within 
an isodensity contour of A = 200, corresponding to a radius r2oo- A 
final cluster catalogue is produced by selecting objects whose total mass 
is equivalent to at least 500 particles of gas and dark matter within this 
radius, and whose centre is not located within a more massive cluster. 
In this way our catalogues (available online from virgo . sussex . ac . uk) 
are complete in mass down to M lim » 1.18 x lO^/i^M©. 

For each cluster the catalogues contain several entries representing 
different cluster properties. Relevant for this study are the bolomet- 
ric luminosity, Lboi = J2 m i Pi A(Ti, ^)/(/xm p ) 2 , the emission-weighted 
temperature, T ew = J2 m i PiA(Tj, Z)Ti/ J2 mi PiA(Tj, Z), the mass- 
weighted temperature T mw = J2 m iTi/ J2 m ii an d the hot gas mass 
fraction, / gas = ^ rrii/M. In these expressions, summations run over 
gas particles with temperatures above 10 5 (1 + z) 2 K (i.e. the hot 
intra-cluster gas that gives rise to the SZ and X-ray emission) within 
a radius r2oo- M is the total mass, i.e., the sum of the mass of all 
particles, including both baryons and dark matter, within the same 
radius. The quantities, mj, Tj, pi and A(Tj, Z) are respectively the mass, 
temperature, density and emissivity (cooling function) of the particles. 

The intrinsic SZ luminosity, defined as Y mt = Fd^fz) 2 , can then 
be computed from the cluster catalogues as, 



where we assume a hydrogen mass fraction of X = 0.76. Since the 
y-parameter is dimensionless we will quote Y mt in units of (/i _1 Mpc) 2 . 
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4. Results: scaling relations at redshift zero 

In Figure 1 we plot the intrinsic SZ-flux against mass (top panel), 
mass-weighted temperature (centre panel) and X-ray luminosity (bot- 
tom panel), for clusters found in simulation boxes at redshift zero. 
All properties are calculated within r2oo- For the plots in this paper we 
adopt the following convention: gray crosses represent clusters from the 
non-radiative simulation, whereas dark gray and light gray crosses illus- 
trate the properties of clusters in the cooling and pre-heating runs. To 
ensure completeness in temperature and luminosity we have trimmed 
the original cluster catalogues in the two bottom plots by selecting 
only clusters with emission-weighted temperatures above 0.35, 0.75 and 
0.55 keV for the non-radiative, cooling and pre-heating simulations, 
respectively. This permits to avoid the introduction of artificial trends 
on the derived fittings for the Ygoo ~~ ^mw an d ^oo ~~ ^boi relations. 

The plots show that Y^qq is tightly correlated with mass and the 
mass-weighted temperature in all runs. The correlation between the 
SZ integrated flux and the X-ray luminosity is also strong, but shows 
significantly more scatter. This is because the X-ray emission is sensi- 
tive to the gas distribution and temperature, whereas Y^oo ^ s essentially 
a function of the total thermal energy of the hot gas mass within r2oo • 
In each plot, the solid lines are power-law best fits to distribution of 
clusters in each simulation. For the Y^oo — -^200 correlation we find the 
following best-fit parameters for each simulation run: 

Non-radiative: Ygjg = 1.2 x 10~ 6 (M 2 oo/Mi 4 ) L7 {h^ 1 Mpc) 2 
Cooling: Ygfc = 8.3 x 10~ 7 (M 2 oo/Mi 4 ) L8 (/T 1 Mpc) 2 

Pre-heating: = 1.0 x 10~ 6 (M 200 /M 14 ) 1 - 7 (h' 1 Mpc) 2 , 

where M14 = 10 14 /i -1 M . The slope from the non-radiative run repro- 
duces well the value 5/3 predicted by the self-similar scaling Eq. (6). 
The cooling simulation presents a steeper slope and lower normalization 
then the other two runs. This implies significant differences in Y^OO' P ar ~ 
ticularly for low-mass systems. For example, at M200 ~ 10 13 /i~ 1 M Q , 
sources in the cooling simulation are about two times less strong then 
their non-radiative counterparts. At 10 14 h~ l M the relative difference 
is still of about 50 percent, but gets much reduced at higher masses. 

The differences between runs can be explained in terms of the changes 
caused by our models of cooling and pre-heating in the hot gas fraction 
and mass-weighted temperature of clusters. The inclusion of cooling in 
both radiative simulations lowers the hot gas mass fraction and tends 
to increases the emission-weighted temperature due to the removal of 
cooled material from the hot gas phase. The increase in temperature 
results from the infall of high-entropy gas which replaces the cooled 
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Figure 1. SZ scaling relations for clusters at redshift zero from the non-radiative 
(gray crosses), cooling (dark gray crosses) and pre-heating (light gray crosses) sim- 
ulations of the ACDM cosmology. From top to bottom we have the scalings of 
the intrinsic F-flux with mass, mass-weighted temperature and X-ray bolometric 
luminosity, respectively. The solid lines are power-law fits to the distributions. 



ws-hcp.tex; 2/02/2008; 10:05; p. 6 



SZ scaling relations in galaxy clusters from hydrodynamical simulations 



7 




Figure 2. Fraction of hot gas within r2oo for clusters in the non-radiative (gray 
crosses) cooling (dark gray crosses) and pre-heating (light gray crosses) simulations. 



material forming at the central regions of clusters. As can be inferred 
from Eq. (7), these are competing effects for the same cluster mass. 
The overall effect on i^oo * s dominated by the decrement of the hot 
gas mass fraction, which shifts clusters downward in the y mt — M 
plane. The effect is less strong in the pre-heating run, because particles 
were heated up by the energy injection mechanism and less amount of 
cooled material is able to form. The decrement of Kjoo * s a ^ so stronger 
in low-mass systems, which explains the steepening of the relations. 

In Figure 2 we plot / gas for all clusters in non-radiative, cooling 
and pre-heating catalogues at redshift zero. We see that clusters in the 
cooling simulation have the lowest hot gas fractions. For low-mass sys- 
tems these are reduced by about 2.5 times when compared with clusters 
in the non-radiative simulation. The fraction of hot gas in clusters of 
the pre-heating run is also substantially reduced, but less than in the 
cooling case. In both radiative simulations, lower mass systems show 
lower amounts of hot gas, which implies a lower Y^oo- ^ n ^ ne coonn g 
simulation, this is because cooling becomes more efficient in low-mass 
systems, where particles have lower temperatures and higher cooling 
rates due to collisional excitation cooling from neutral elements. In the 
pre-heating simulation, low-mass systems may also reduce their hot 
gas fraction in this way. This process is less efficient then in the cooling 
case because the energy injection rises the temperature of the particles 
and reduces the contribution from collisional excitation cooling. Pre- 
heating can also reduce the amount of hot gas in clusters in a different 
way. Instead of converting hot gas into cooled material, pre-heating has 
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the effect of heating up and expelling gas from the central regions of 
clusters. In low-mass systems this may cause some of the hot gas to 
be expelled from the haloes, which in turn implies lower / gas at lower 
masses. We find that the amount of baryonic gas which is converted into 
cooled material in clusters of the cooling and pre-heating simulation can 
be as high as 50 and 30 per cent, respectively. 

The correlation between Y^oo an d the mass-weighted temperature is 
also easy to understand. For the same T mw , Eq. (7) tells us that the Y- 
luminosity is simply proportional to the hot gas mass of the cluster. As 
a consequence, for a given temperature, clusters of the cooling and pre- 
heating simulations present lower Y^oo- The effect is again stronger for 
low-mass systems and for clusters in the cooling simulation. Clusters 
in the radiative runs also present higher temperatures for the same 
mass. The combination of these two effects causes clusters to shift in 
the direction of high temperatures and lower Y-fluxes in the Y^qq — T mw 
plane. The best-fit parameters found for the Y^qq — T mw relation are: 
Non-radiative: = 5.8 x lO" 7 (T mw /keV) 2 - 6 (h' 1 Mpc) 2 

Cooling: = 1.8 x 10~ 7 (T mw /keV) 3 - (/i^ 1 Mpc) 2 

Pre-heating: = 3.4 x 10~ 7 (T mw /keV) 2 - 8 (h^ 1 Mpc) 2 

The slopes of the relations derived from the radiative simulations are 
steeper than those found in the non-radiative run. This is due to the 
higher temperatures and lower gas masses of clusters in these simula- 
tions. Again the non-radiative simulation is in good agreement with 
the slope 5/2 predicted by the self-similar scaling. 

In the top panels of Figure 3 we plot Y^qq against the mass- weighted 
(left panel) and X-ray emission-weighted (right panel) temperatures 
for clusters in the non-radiative and cooling simulations. These show 
that the integrated SZ signal presents a much tighter correlation with 
the mass-weighted temperature than with the emission-weighted tem- 
perature. This result confirms expectations, since Tew is weighted by 
density and temperature and is therefore sensitive to substructures and 
clumping of the gas inside the cluster. The bottom panels of Figure 3 
show the correlation between mass and temperature for the same set 
of clusters displayed in the top panels. Again, the correlation between 
M200 and T mw presents significantly less scatter than the correlation of 
-^200 with T ew . This is evidence of the potential of the SZ effect to study 
the properties of clusters. Future high-resolution SZ observations may 
provide important constraints on cluster evolutionary theories. Good 
calibrations of the observed Y—T and Y—M relations may also provide 
alternative robust estimates of temperature and mass in clusters. 

We end the analysis by reporting on the Y^qq— L\ m \ correlation. In the 
radiative runs, the high entropy gas which replaces the cooled material 
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Figure 3. The top panels show the Y^'oo against the mass-weighted (left panel) and 
X-ray emission-weighted (right panel) temperatures (respectively T mw and T ew ) for 
clusters in the non-radiative (gray crosses) and cooling (dark gray crosses) runs. The 
bottom panels show mass-temperature correlations for the same set of clusters. 



which forms in these simulations is less dense and hotter then the gas 
in the non-radiative simulation. These variations in temperature and 
density have opposite effects on the X-ray luminosity, which is a slowly- 
varying function of temperature, but it is proportional to the square of 
the gas density, see Eq. (3). As a result, the change in density wins the 
competition and the luminosity decreases with cooling. The inclusion 
of pre-heating reduces the amount of cooled material and increases the 
gas temperature. This implies higher luminosities in the pre-heating 
run then in the cooling simulation, but lower luminosities than in the 
non-radiative run. The combination of these effects with the decrease of 
Y200 explains the substantial shifts observed in the Y 2 m o ~~ ^boi relations 
of Figure 1 in the pre-heating and cooling simulations. The best-fit 
parameters for the distribution of clusters in the — Lbol plane are: 
Non-radiative: Kjoo = 6 - 3 x 10 ~ 12 ( L boi/U ) 13 (h' 1 Mpc) 2 
Cooling: Y$& = 1.6 x 1CT 9 (L hol /L m f 9 {h' 1 Mpc) 2 

Pre-heating: = 3.2 x 1CT 10 (L bol /L 40 ) L0 (hr 1 Mpc) 2 , 

where L40 = 10 40 /i _2 erg s . These show considerable differences of 
slope and normalization. The slope in the non-radiative run is close 
to 1.25, which is the value predicted by the self-similar scaling. The 
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radiative runs show similar slopes, which are less steep than in the non- 
radiative case. As in the case of the Kjqq — T cw relation, the Kjoo ~~ A>oi 
relation shows significant scatter due to the sensitivity of the X-ray 
emission to the density and temperature distribution of the hot gas. 

5. Conclusions 

We have studied the effect of our models of cooling and pre-heating on 
the SZ cluster population by identifying clusters in the simulation boxes 
and computing their characteristic properties. We correlated the SZ 
luminosity with other cluster properties and derived scaling relations at 
redshift zero. The non-radiative simulation reproduces well the scalings 
predicted by the self-similar model, whereas the inclusion of cooling 
and pre-heating generally changes the slope and normalization of the 
scaling relations for those runs. The integrated Y-signal is found to be 
tightly correlated with the mass-weighted temperature and total mass 
of the cluster. 

The author wishes to thank A. Liddle, P. Thomas, S. Kay and 
O. Muanwong for many fruitful discussions and acknowledges the use 
of the computer facilities of the Astronomy Unit at Sussex, UK, and 
CALMIP at Toulouse, France. The simulations used were carried out on 
a Cray-T3E at EPCC as part of the VIRGO Consortium collaboration. 
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